* Wealth index calculation. Claudia Marchena, 01/18/2017. Ghana MIS 2016. cd 'C:\users\21762\Root of C\hnp2a\Ghana MIS 2016'. include file='C:\users\21762\Root of C\hnp2a\Ghana MIS 2016\wealth_assets.sps'. execute. set tvars=both. * Macro definitions to perform some recoding automatically. * Note: Do not use an asterisk in front of a macro name to try to comment it out - it will run any way because of the way the macros are expanded. * Delete the line instead or change the spelling of the macro in the comment. * Do not use the names of macros in comments - they will also be expanded too, and likely will fail. *{Create binary variables based on condition and give label }. define dichotomize( !positional !enclose('(',')') / !positional !tokens(1) / !positional !tokens(1) ) compute !2=0. if (!1) !2=1. variable labels !2 !3. value labels !2 0 'No' 1 'Yes'. formats !2 (f1.0). !enddefine. *(Separate shared toilets from non-shared toilets). define shared_fac( !positional !tokens(1) / !positional !tokens(1) / !positional !tokens(1) ). compute !2 = 0. do if (QH106=1). + if (!1=1) !2 = 1. + compute !1 = 0. end if. variable labels !2 !3. value labels !2 0 'No' 1 'Yes'. formats !2 (f1.0). !enddefine. *{Reset missing values to 'does not have', change 2 code to 0}. define no2zero( !positional !tokens(1) ) if (missing(!1) | !1<>1) !1=0. value labels !1 0 'No' 1 'Yes'. !enddefine. *{Construct Variables}. * Source of drinking water. dichotomize (QH101=11) QH101_11 "Source of drinking water: Piped into dwelling". dichotomize (QH101=12) QH101_12 "Source of drinking water: Piped to yard/plot". dichotomize (QH101=13) QH101_13 "Source of drinking water: Piped to neighbor". dichotomize (QH101=14) QH101_14 "Source of drinking water: Public tap/standpipe". dichotomize (QH101=21) QH101_21 "Source of drinking water: Tube well or borehole". dichotomize (QH101=31) QH101_31 "Source of drinking water: Protected well". dichotomize (QH101=32) QH101_32 "Source of drinking water: Unprotected well". dichotomize (QH101=41) QH101_41 "Source of drinking water: Protected spring". dichotomize (QH101=42) QH101_42 "Source of drinking water: Unprotected spring". dichotomize (QH101=51) QH101_51 "Source of drinking water: Rainwater". dichotomize (QH101=61) QH101_61 "Source of drinking water: Tanker truck". dichotomize (QH101=71) QH101_71 "Source of drinking water: Cart with small tank". dichotomize (QH101=81) QH101_81 "Source of drinking water: Surface water (river/dam/lake/pond/stream/canal/irrigation channel)". dichotomize (QH101=91) QH101_91 "Source of drinking water: Bottled water". dichotomize (QH101=92) QH101_92 "Source of drinking water: Sachet water". dichotomize (QH101=96) QH101_96 "Source of drinking water: Other". * Type of toilet facility. dichotomize (QH105=11) QH105_11 "Type of toilet facility: Flush to piped sewer system". dichotomize (QH105=12) QH105_12 "Type of toilet facility: Flush to septic tank". dichotomize (QH105=13) QH105_13 "Type of toilet facility: Flush to pit latrine". dichotomize (QH105=14) QH105_14 "Type of toilet facility: Flush to somewhere else". dichotomize (QH105=15) QH105_15 "Type of toilet facility: Flush, don't know where". dichotomize (QH105=21) QH105_21 "Type of toilet facility: Ventilated improved pit latrine". dichotomize (QH105=22) QH105_22 "Type of toilet facility: Pit latrine with slab". dichotomize (QH105=23) QH105_23 "Type of toilet facility: Pit latrine without slab/open pit". dichotomize (QH105=31) QH105_31 "Type of toilet facility: Composting toilet". dichotomize (QH105=41) QH105_41 "Type of toilet facility: Bucket toilet". dichotomize (QH105=51) QH105_51 "Type of toilet facility: Hanging toilet/hanging latrine". dichotomize (QH105=61) QH105_61 "Type of toilet facility: No facility/bush/field". dichotomize (QH105=96) QH105_96 "Type of toilet facility: Other". shared_fac QH105_11 QH105_11_sh "Type of toilet facility: Flush to piped sewer system - shared". shared_fac QH105_12 QH105_12_sh "Type of toilet facility: Flush to septic tank - shared". shared_fac QH105_13 QH105_13_sh "Type of toilet facility: Flush to pit latrine - shared". shared_fac QH105_14 QH105_14_sh "Type of toilet facility: Flush to somewhere else - shared". shared_fac QH105_15 QH105_15_sh "Type of toilet facility: Flush, don't know where - shared". shared_fac QH105_21 QH105_21_sh "Type of toilet facility: Ventilated improved pit latrine - shared". shared_fac QH105_22 QH105_22_sh "Type of toilet facility: Pit latrine with slab - shared". shared_fac QH105_23 QH105_23_sh "Type of toilet facility: Pit latrine without slab/open pit - shared". shared_fac QH105_31 QH105_31_sh "Type of toilet facility: Composting toilet - shared". shared_fac QH105_41 QH105_41_sh "Type of toilet facility: Bucket toilet - shared". shared_fac QH105_51 QH105_51_sh "Type of toilet facility: Hanging toilet/hanging latrine - shared". shared_fac QH105_96 QH105_96_sh "Type of toilet facility: Other - shared". * Type of cooking fuel. dichotomize (QH108=1) QH108_1 "Type of cooking fuel: Electricity". dichotomize (QH108=2) QH108_2 "Type of cooking fuel: LPG". dichotomize (QH108=3) QH108_3 "Type of cooking fuel: Natural gas". dichotomize (QH108=4) QH108_4 "Type of cooking fuel: Biogas". dichotomize (QH108=5) QH108_5 "Type of cooking fuel: Kerosene". dichotomize (QH108=6) QH108_6 "Type of cooking fuel: Coal, lignite". dichotomize (QH108=7) QH108_7 "Type of cooking fuel: Charcoal". dichotomize (QH108=8) QH108_8 "Type of cooking fuel: Wood". dichotomize (QH108=9) QH108_9 "Type of cooking fuel: Straw/shrubs/grass". dichotomize (QH108=10) QH108_10 "Type of cooking fuel: Agricultural crop residue". dichotomize (QH108=11) QH108_11 "Type of cooking fuel: Animal dung". dichotomize (QH108=95) QH108_95 "Type of cooking fuel: No food cooked in household". dichotomize (QH108=96) QH108_96 "Type of cooking fuel: Other". * Cows/bulls. if (missing(QH111A) | QH110 <> 1) QH111A=0. missing values QH111A (99). compute QH111A_0 = (QH111A = 0). compute QH111A_1 = (QH111A >= 1 & QH111A <= 4). compute QH111A_2 = (QH111A >= 5 & QH111A <= 9). compute QH111A_3 = (QH111A >= 10 & QH111A <= 98). variable labels QH111A_0 'Cows/bulls: None' /QH111A_1 'Cows/bulls: 1-4' /QH111A_2 'Cows/bulls: 5-9' /QH111A_3 'Cows/bulls: 10+'. * Other cattle. if (missing(QH111B) | QH110 <> 1) QH111B=0. missing values QH111B (99). compute QH111B_0 = (QH111B = 0). compute QH111B_1 = (QH111B >= 1 & QH111B <= 4). compute QH111B_2 = (QH111B >= 5 & QH111B <= 9). compute QH111B_3 = (QH111B >= 10 & QH111B <= 98). variable labels QH111B_0 'Other cattle: None' /QH111B_1 'Other cattle: 1-4' /QH111B_2 'Other cattle: 5-9' /QH111B_3 'Other cattle: 10+'. * Horses/donkeys/mules. if (missing(QH111C) | QH110 <> 1) QH111C=0. missing values QH111C (99). compute QH111C_0 = (QH111C = 0). *compute QH111C_1 = (QH111C >= 1 & QH111C <= 4). *compute QH111C_2 = (QH111C >= 5 & QH111C <= 9). *compute QH111C_3 = (QH111C >= 10 & QH111C <= 98). *variable labels QH111C_0 'Horses/donkeys/mules: None' /QH111C_1 'Horses/donkeys/mules: 1-4' /QH111C_2 'Horses/donkeys/mules: 5-9' /QH111C_3 'Horses/donkeys/mules: 10+'. compute QH111C_1 = (QH111C > 0 and qh111c <=95). variable labels QH111C_0 'Horses/donkeys/mules: None' /QH111C_1 'Horses/donkeys/mules: 1+'. * Goats. if (missing(QH111D) | QH110 <> 1) QH111D=0. missing values QH111D (99). compute QH111D_0 = (QH111D = 0). compute QH111D_1 = (QH111D >= 1 & QH111D <= 4). compute QH111D_2 = (QH111D >= 5 & QH111D <= 9). compute QH111D_3 = (QH111D >= 10 & QH111D <= 98). variable labels QH111D_0 'Goats: None' /QH111D_1 'Goats: 1-4' /QH111D_2 'Goats: 5-9' /QH111D_3 'Goats: 10+'. * Sheep. if (missing(QH111E) | QH110 <> 1) QH111E=0. missing values QH111E (99). compute QH111E_0 = (QH111E = 0). compute QH111E_1 = (QH111E >= 1 & QH111E <= 4). compute QH111E_2 = (QH111E >= 5 & QH111E <= 9). compute QH111E_3 = (QH111E >= 10 & QH111E <= 98). variable labels QH111E_0 'Sheep: None' /QH111E_1 'Sheep: 1-4' /QH111E_2 'Sheep: 5-9' /QH111E_3 'Sheep: 10+'. * Chickens or other poultry. if (missing(QH111F) | QH110 <> 1) QH111F=0. missing values QH111F (99). compute QH111F_0 = (QH111F = 0). compute QH111F_1 = (QH111F >= 1 & QH111F <= 9). compute QH111F_2 = (QH111F >= 10 & QH111F <= 29). compute QH111F_3 = (QH111F >= 30 & QH111F <= 98). variable labels QH111F_0 'Chickens or other poultry: None' /QH111F_1 'Chickens or other poultry: 1-4' /QH111F_2 'Chickens or other poultry: 5-9' /QH111F_3 'Chickens or other poultry: 10+'. * Pigs. if (missing(QH111G) | QH110 <> 1) QH111G=0. missing values QH111G (99). compute QH111G_0 = (QH111G = 0). compute QH111G_1 = (QH111G >= 1 & QH111G <= 4). compute QH111G_2 = (QH111G >= 5 & QH111G <= 9). compute QH111G_3 = (QH111G >= 10 & QH111G <= 98). variable labels QH111G_0 'Pigs: None' /QH111G_1 'Pigs: 1-4' /QH111G_2 'Pigs: 5-9' /QH111G_3 'Pigs: 10+'. * Rabbits. if (missing(QH111H) | QH110 <> 1) QH111H=0. missing values QH111H (99). compute QH111H_0 = (QH111H = 0). *compute QH111H_1 = (QH111H >= 1 & QH111H <= 9). *compute QH111H_2 = (QH111H >= 10 & QH111H <= 29). *compute QH111H_3 = (QH111H >= 20 & QH111H <= 98). *variable labels QH111H_0 'Rabbits: None' /QH111H_1 'Rabbits: 1-4' /QH111H_2 'Rabbits: 5-9' /QH111H_3 'Rabbits: 10+'. compute QH111H_1 = (QH111H > 0 and qh111H <=95). variable labels QH111H_0 'Rabbits: None' /QH111H_1 'Rabbits: 1+'. * Grasscutter. if (missing(QH111I) | QH110 <> 1) QH111I=0. missing values QH111I (99). compute QH111I_0 = (QH111I = 0). *compute QH111I_1 = (QH111I >= 1 & QH111I <= 4). *compute QH111I_2 = (QH111I >= 10 & QH111I <= 29). *compute QH111I_3 = (QH111I >= 20 & QH111I <= 98). *variable labels QH111I_0 'Grasscutter: None' /QH111I_1 'Grasscutter: 1-4' /QH111I_2 'Grasscutter: 5-9' /QH111I_3 'Grasscutter: 10+'. compute QH111I_1 = (QH111I > 0 and qh111I <=95). variable labels QH111I_0 'Grasscutter: None' /QH111I_1 'Grasscutter: 1+'. * Hectares for agricultural land (number). compute landarea = 0. if (QH113U = 1) landarea = QH113N. /* hectares */ if (QH113U = 2) landarea = QH113N*(0.404686). /* 1 Acre = 0.404686 hectares */ if (QH113U = 3) landarea = QH113N*(0.101171). /* 1 plot = 0.101171 hectares */ if (QH113U = 9 and QH113N = 95.0) landarea = 95 * 0.404686. *if (QH113U >= 8) landarea = 0. if (missing(QH113N) | QH113N >= 99.8) landarea=999999. if (missing(QH112 ) | QH112 <> 1) landarea=0. missing values landarea (999999). * Electricity. no2zero QH114A. * Radio. no2zero QH114B. * Television. no2zero QH114C. * Telephone (non-mobile). no2zero QH114D. * Computer. no2zero QH114E. * Refrigerator. no2zero QH114F. * Freezer. no2zero QH114G. * Generator. no2zero QH114H. * Washing machine. no2zero QH114I. * Camera. no2zero QH114J. * Video/DVD/VCD. no2zero QH114K. * Sewing machine. no2zero QH114L. * Bed. no2zero QH114M. * Table. no2zero QH114N. * Chair. no2zero QH114O. * Cabinet. no2zero QH114P. * Wrist watch. no2zero QH115A. * Mobile telephone. no2zero QH115B. * Bicycle. no2zero QH115C. * Motorcycle or scooter. no2zero QH115D. * Animal-drawn cart. no2zero QH115E. * Car or Truck. no2zero QH115F. * Boat with a motor. no2zero QH115G. * Boat without a motor. no2zero QH115H. * Bank account. no2zero QH116. *{Members per sleeping room}. if (hhusual=0) hhusual=hhslept. if (QH109>0) memsleep=trunc(hhusual/QH109). if (QH109=0) memsleep=hhusual. if (missing(QH109) or QH109>=99 or memsleep>=98) memsleep=99. variable labels memsleep 'Number of members per sleeping room'. value labels memsleep 0 'Less than 1 per room'. formats memsleep (f2.0). missing values memsleep (99). * Compute urban and rural variables coded (1/0) for filters later. compute urban=(QHTYPE = 1). compute rural=(QHTYPE = 2). variable labels urban 'Urban' / rural 'Rural'. value labels urban 1 'Urban' / rural 1 'Rural'. formats urban rural (f1.0). execute. * Check on indicator variable creation. frequencies variables=QHTYPE to domestic. * Toilet facility by shared/not shared. crosstabs /tables=QH105 by QH106. * Land area by units - if there are separate units - need to convert them to one unit. crosstabs /tables=QH113N by QH113U. frequencies variables= QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_92 QH101_96 QH105_11 QH105_12 QH105_13 QH105_14 QH105_15 QH105_21 QH105_22 QH105_23 QH105_31 QH105_41 QH105_51 QH105_61 QH105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_15_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_31_sh QH105_41_sh QH105_51_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_5 QH108_6 QH108_7 QH108_8 QH108_9 QH108_10 QH108_11 QH108_95 QH108_96 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH114N QH114O QH114P QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH115H QH116 DOMESTIC HOUSE LAND memsleep QH111A_0 QH111A_1 QH111A_2 QH111A_3 QH111B_0 QH111B_1 QH111B_2 QH111B_3 QH111C_0 QH111C_1 QH111D_0 QH111D_1 QH111D_2 QH111D_3 QH111E_0 QH111E_1 QH111E_2 QH111E_3 QH111F_0 QH111F_1 QH111F_2 QH111F_3 QH111G_0 QH111G_1 QH111G_2 QH111G_3 QH111H_0 QH111H_1 QH111I_0 QH111I_1 landarea . * Turn off weights before all factor analysis. weight off. * Name the dataset window for the HH data for use later. dataset name assets. save outfile='assets.sav'. *** Test Factor Analysis. filter off. execute. * removed due to no cases: qh101_96 qh105_14 qh105_15 qh105_31 qh108_4 qh108_11 qh108_96 qh111c_2. * recode few cases. if (qh105_41 eq 1 or qh105_51 eq 1) qh105_96=1. if (qh105_31_sh eq 1 or qh105_41_sh eq 1 or qh105_51_sh eq 1) qh105_96_sh=1. if (qh108_3 eq 1) qh108_2=1. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_21 QH105_22 QH105_23 QH105_61 qh105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_15_sh QH105_21_sh QH105_22_sh QH105_23_sh qh105_96_sh QH108_1 QH108_2 QH108_5 QH108_6 QH108_7 QH108_8 QH108_9 QH108_10 QH108_95 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH114N QH114O QH114P QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH115H QH116 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_2 QH111B_3 QH111C_1 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_0 QH111G_1 QH111G_2 QH111G_3 QH111H_0 QH111H_1 QH111I_0 QH111I_1 landarea /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_21 QH105_22 QH105_23 QH105_61 qh105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_15_sh QH105_21_sh QH105_22_sh QH105_23_sh qh105_96_sh QH108_1 QH108_2 QH108_5 QH108_6 QH108_7 QH108_8 QH108_9 QH108_10 QH108_95 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH114N QH114O QH114P QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH115H QH116 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_2 QH111B_3 QH111C_1 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_0 QH111G_1 QH111G_2 QH111G_3 QH111H_0 QH111H_1 QH111I_0 QH111I_1 landarea /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /method=correlation. *** Common Factor analysis ***. ** Removed area-specific variables **. ** Agricultural animal variables excluded ** Land area excluded ** Any others ?. *** Common Factor Analysis. filter off. execute. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_21 QH105_22 QH105_23 QH105_61 qh105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_15_sh QH105_21_sh QH105_22_sh QH105_23_sh qh105_96_sh QH108_1 QH108_2 QH108_5 QH108_6 QH108_7 QH108_8 QH108_9 QH108_10 QH108_95 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH114N QH114O QH114P QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH115H QH116 memsleep /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_21 QH105_22 QH105_23 QH105_61 qh105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_15_sh QH105_21_sh QH105_22_sh QH105_23_sh qh105_96_sh QH108_1 QH108_2 QH108_5 QH108_6 QH108_7 QH108_8 QH108_9 QH108_10 QH108_95 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH114N QH114O QH114P QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH115H QH116 memsleep /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /save reg(ALL com) /method=correlation. *** Urban Factor Analysis. filter off. filter by urban. execute. * removed due to no cases: qh111i_3. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_21 QH105_22 QH105_23 QH105_61 qh105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_15_sh QH105_21_sh QH105_22_sh QH105_23_sh qh105_96_sh QH108_1 QH108_2 QH108_5 QH108_6 QH108_7 QH108_8 QH108_9 QH108_10 QH108_95 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH114N QH114O QH114P QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH115H QH116 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_2 QH111B_3 QH111C_1 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_0 QH111G_1 QH111G_2 QH111G_3 QH111H_0 QH111H_1 QH111I_0 QH111I_1 landarea /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_21 QH105_22 QH105_23 QH105_61 qh105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_15_sh QH105_21_sh QH105_22_sh QH105_23_sh qh105_96_sh QH108_1 QH108_2 QH108_5 QH108_6 QH108_7 QH108_8 QH108_9 QH108_10 QH108_95 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH114N QH114O QH114P QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH115H QH116 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_2 QH111B_3 QH111C_1 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_0 QH111G_1 QH111G_2 QH111G_3 QH111H_0 QH111H_1 QH111I_0 QH111I_1 landarea /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /save reg(ALL urb) /method=correlation. *** Rural Factor Analysis. filter off. filter by rural. execute. * removed due to no cases: qh105_15_sh. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_21 QH105_22 QH105_23 QH105_61 qh105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_21_sh QH105_22_sh QH105_23_sh qh105_96_sh QH108_1 QH108_2 QH108_5 QH108_6 QH108_7 QH108_8 QH108_9 QH108_10 QH108_95 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH114N QH114O QH114P QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH115H QH116 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_2 QH111B_3 QH111C_1 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_0 QH111G_1 QH111G_2 QH111G_3 QH111H_0 QH111H_1 QH111I_0 QH111I_1 landarea /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_92 QH105_11 QH105_12 QH105_13 QH105_21 QH105_22 QH105_23 QH105_61 qh105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_21_sh QH105_22_sh QH105_23_sh qh105_96_sh QH108_1 QH108_2 QH108_5 QH108_6 QH108_7 QH108_8 QH108_9 QH108_10 QH108_95 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH114N QH114O QH114P QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH115H QH116 LAND memsleep QH111A_1 QH111A_2 QH111A_3 QH111B_1 QH111B_2 QH111B_3 QH111C_1 QH111D_1 QH111D_2 QH111D_3 QH111E_1 QH111E_2 QH111E_3 QH111F_1 QH111F_2 QH111F_3 QH111G_0 QH111G_1 QH111G_2 QH111G_3 QH111H_0 QH111H_1 QH111I_0 QH111I_1 landarea /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /save reg(ALL rur) /method=correlation. * Label the created score variables. rename variables (com1 urb1 rur1=comscore urbscore rurscore). variable labels comscore 'Common wealth score' /urbscore 'Urban wealth score' /rurscore 'Rural wealth score'. * Add a variable used for linking later. filter off. string rowtype_ (A8). compute rowtype_ = 'EST'. * Calculate regressions. ** Area=urban. filter off. filter by urban. execute. * Declare a dataset to be written to in the regression. dataset declare urbancorv. * Run regression of comscore with urbscore. regression /missing listwise /statistics coeff outs R anova /criteria=pin(.05) pout(.10) /noorigin /dependent comscore /method=enter urbscore /outfile=corv(urbancorv). * Activate file of output from regression. dataset activate urbancorv. * Drop all rows of output except the coefficients. select if (rowtype_ = 'EST'). execute. * Delete unnecessary variables before merging. delete variables DEPVAR_ VARNAME_. * Rename variables containing the constant and the coefficient. rename variables CONST_=urbconst urbscore=urbcoeff. * Re-activate the main household data. dataset activate assets. * merge the coefficients. match files /file = * /table = urbancorv /by ROWTYPE_. execute. * Calculate regressions. ** Area=rural. filter off. filter by rural. execute. * Declare a dataset to be written to in the regression. dataset declare ruralcorv. * Run regression of comscore with rurscore. regression /missing listwise /statistics coeff outs R anova /criteria=pin(.05) pout(.10) /noorigin /dependent comscore /method=enter rurscore /outfile=corv(ruralcorv). * Activate file of output from regression. dataset activate ruralcorv. * Drop all rows of output except the coefficients. select if (rowtype_ = 'EST'). execute. * Delete unnecessary variables before merging. delete variables DEPVAR_ VARNAME_. * Rename variables containing the constant and the coefficient. rename variables CONST_=rurconst rurscore=rurcoeff. * Re-activate the main household data. dataset activate assets. * merge the coefficients. match files /file = * /table = ruralcorv /by ROWTYPE_. execute. dataset close urbancorv. dataset close ruralcorv. dataset activate assets. filter off. *** Calculate combined wealth score from Urban and Rural Scores. * Use coefficients from urban and rural regressions above. compute combscor=0. print formats combscor (f11.5). write formats combscor (f11.5). ** Urban. if (urban = 1) combscor=urbconst+urbcoeff*urbscore. ** Rural. if (rural = 1) combscor=rurconst+rurcoeff*rurscore. variable labels combscor 'Combined national wealth score'. execute. compute hhwt = QHWEIGHT/1000000. variable labels hhwt 'HH weights'. weight by hhwt. filter off. execute. frequencies variables=combscor comscore /format=notable /ntiles=5 /statistics=stddev mean /histogram normal /order=analysis. filter off. filter by urban. execute. frequencies variables=combscor urbscore /format=notable /ntiles=5 /statistics=stddev mean /histogram normal /order=analysis. filter off. filter by rural. execute. frequencies variables=combscor rurscore /format=notable /ntiles=5 /statistics=stddev mean /histogram normal /order=analysis. filter off. execute. *Calculate quintiles and scores for data file. compute hhmemwt=QHWEIGHT*hhusual/1000000. weight by hhmemwt. variable labels hhmemwt 'HH members weighting for index'. filter off. filter by urban. execute. rank variables=urbscore (A) /rank /ntiles (5) /print=yes /ties=mean. filter off. filter by rural. execute. rank variables=rurscore (A) /rank /ntiles (5) /print=yes /ties=mean. filter off. execute. rank variables=combscor (A) /rank /ntiles (5) /print=yes /ties=mean. variable labels ncombsco 'Combined wealth index' /nurbscor 'Urban wealth index' /nrurscor 'Rural wealth index'. value labels ncombsco nurbscor nrurscor 1 'Lowest' 2 'Second' 3 'Middle' 4 'Fourth' 5 'Highest'. *** Check on quintiles. frequencies variables=ncombsco nurbscor nrurscor. weight by hhwt. ctables /table ( QH101_11+QH101_12+QH101_13+QH101_14+QH101_21+QH101_31+QH101_32+QH101_41+QH101_42+QH101_51+QH101_61+QH101_71+QH101_81+QH101_91+QH101_92+QH101_96+ QH105_11+QH105_12+QH105_13+QH105_14+QH105_15+QH105_21+QH105_22+QH105_23+QH105_31+QH105_41+QH105_51+QH105_61+QH105_96+ QH105_11_sh+QH105_12_sh+QH105_13_sh+QH105_14_sh+QH105_15_sh+QH105_21_sh+QH105_22_sh+QH105_23_sh+QH105_31_sh+QH105_41_sh+QH105_51_sh+QH105_96_sh+ QH108_1+QH108_2+QH108_3+QH108_4+QH108_5+QH108_6+QH108_7+QH108_8+QH108_9+QH108_10+QH108_11+QH108_95+QH108_96+ QH114A+QH114B+QH114C+QH114D+QH114E+QH114F+QH114G+QH114H+QH114I+QH114J+QH114K+QH114L+QH114M+QH114N+QH114O+QH114P+QH115A+QH115B+QH115C+QH115D+QH115E+QH115F+ QH115G+QH115H+QH116+DOMESTIC+HOUSE+LAND+memsleep+ QH111A_0+QH111A_1+QH111A_2+QH111A_3+QH111B_0+QH111B_1+QH111B_2+QH111B_3+ QH111C_0+QH111C_1+QH111C_2+QH111C_3+QH111D_0+QH111D_1+QH111D_2+QH111D_3+ QH111E_0+QH111E_1+QH111E_2+QH111E_3+QH111F_0+QH111F_1+QH111F_2+QH111F_3+ QH111G_0+QH111G_1+QH111G_2+QH111G_3+QH111H_0+QH111H_1+QH111H_2+QH111H_3+ QH111I_0+QH111I_1+QH111I_2+QH111I_3+landarea )[S][Mean F8.3] by (ncombsco+nurbscor+nrurscor)[C] /slabels visible=no. * Alternative instead of ctables. *means tables= QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_31 QH101_32 QH101_41 QH101_42 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_92 QH101_96 QH105_11 QH105_12 QH105_13 QH105_14 QH105_15 QH105_21 QH105_22 QH105_23 QH105_31 QH105_41 QH105_51 QH105_61 QH105_96 QH105_11_sh QH105_12_sh QH105_13_sh QH105_14_sh QH105_15_sh QH105_21_sh QH105_22_sh QH105_23_sh QH105_31_sh QH105_41_sh QH105_51_sh QH105_96_sh QH108_1 QH108_2 QH108_3 QH108_4 QH108_5 QH108_6 QH108_7 QH108_8 QH108_9 QH108_10 QH108_11 QH108_95 QH108_96 QH114A QH114B QH114C QH114D QH114E QH114F QH114G QH114H QH114I QH114J QH114K QH114L QH114M QH114N QH114O QH114P QH115A QH115B QH115C QH115D QH115E QH115F QH115G QH115H QH116 DOMESTIC HOUSE LAND memsleep QH111A_0 QH111A_1 QH111A_2 QH111A_3 QH111B_0 QH111B_1 QH111B_2 QH111B_3 QH111C_0 QH111C_1 QH111C_2 QH111C_3 QH111D_0 QH111D_1 QH111D_2 QH111D_3 QH111E_0 QH111E_1 QH111E_2 QH111E_3 QH111F_0 QH111F_1 QH111F_2 QH111F_3 QH111G_0 QH111G_1 QH111G_2 QH111G_3 QH111H_0 QH111H_1 QH111H_2 QH111H_3 QH111I_0 QH111I_1 QH111I_2 QH111I_3 landarea by ncombsco nurbscor nrurscor /cells mean count stddev. graph /histogram(normal)=combscor /title= 'Distribution of Households by Wealth Scores'. frequencies variables=combscor /format=notable /ntiles=5 /statistics=stddev minimum maximum semean mean median mode skewness seskew kurtosis sekurt /order=analysis. weight off. filter off. write formats combscor urbscore rurscore (f11.5). * Save final dataset of assets. save outfile='assets.sav'. *** Write out scores file. write outfile='scores.dat' table /QHCLUST QHNUMBER combscor ncombsco urbscore nurbscor rurscore nrurscor. execute.